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Abstract 

A new general .fT-network reliability factorization theorem is proved. Beside 
its theoretical mathematical importance the theorem shows how to do parallel 
processing in exact network reliability calculations in order to reduce the 
processing time of this NP-hard problem. It also shows a factorization of the 
Random Cluster Model of Statistical Mechanics in the non cluster limit. 

Keywords: Network Reliability, Graph Theory, Factorization, Parallel 
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1. Introduction 

The i^-network reliability factorization theorem [Mo] and the reduction 
transformations (series-parallel, polygon-to-chain [Wo] and delta-star [Ga]) 
are the key stone of the known factoring algorithms [SC] for the network 
reliability exact calculation. Beside these results and the well known factor- 
ization through an articulation point, no other general factorization theorem 
is known in exact i^-network reliability calculation. This paper gives a new 
general factorization theorem solving the following problem: 

Problem: Given a decomposition of a stochastic graph G by subgraphs 
G\ and G 2 only sharing nodes, express the reliability of G in terms of the 
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reliabilities of the graphs resulting from G\ and G 2 identifying the common 
nodes shared by them in all possible ways. 

Exact calculation of fT-network reliability is a NP-hard problem [Co] so it 
is necessary to know how to do parallel processing in order to reduce the pro- 
cessing time. This theorem has the following immediate parallel processing 
application: Suppose we are willing to calculate the exact network reliability 
of a stochastic graph G = (V, E, (p e ) e es) with m edges and perfect nodes. A 
priori the processing time involve 2 m steps. Cutting the graph G through a 
set of nodes obtaining G\ and G2 with approximately m/2 edges each, the 
solution of the above problem allow us to calculate the reliability of G with 
2 m / 2 steps in each processing line. 

As another application, it is shown in the appendix a non trivial fac- 
torization of the Random Cluster Model of statistical mechanics in the non 
cluster limit. 

2. Preliminaries 

2.1. K- Reliability 

The mathematical model of a Network whose nodes are perfect and its 
edges can fail is a stochastic graph [Co]; i.e. an undirected graph with asso- 
ciated Bernoulli variables to its edges. 

Definition 2.1. An undirected graph G is a pair (V, E) such that V is finite 
set whose elements will be called nodes and E is a subset of {{a, b} / a,b E 
V} x N whose elements will be called edges such that for each pair of distinct 
edges ({oti, 61}, rii), ({02, b 2 }, n 2 ) E E we have that n\ 7^ n 2 . 

Definition 2.2. A stochastic graph G is a tern (V, E, $) such that (V, E) is 
a graph and $ : E —¥ Ber is a function which associates a Bernoulli variable 
to each edge in such a way that these variables are independent. 

Each Bernoulli variable is characterized by a parameter p in the [0, 1] 
closed interval and we can write a stochastic graph as (G, {p e }eeE) where G 
is an undirected graph and p e is the parameter of the variable $(e). Nodes 
and edges of G will be denoted by V{G) and E(G) respectively. 

Definition 2.3. A state E of the graph G = (V, E) is a function 8 : E — >■ 
{0, 1}. An edge e will be called operative if £(e) = 1 and will be called 
non-operative otherwise. 
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Consider a subset K of V(G). A state £ of the graph G will be called 
a f^-PathSet (or inoperative) if K is contained in the set of nodes of any 
of the edge-connected components of the graph resulting from removing the 
non-operative edges of G. Otherwise the state will be called a If-CutSet. 

Definition 2.4. The K- Reliability of a stochastic graph G is 

R K {G) = P(£ is a K — PathS et) 

Because of the independence of the bernoulli variables associated to the 
edges, we can calculate the i^-Reliability in the following way: 

p{£)= n Pi^^-Pi) 1 -"^ 

<%eE(G) 

R K {G) = p ( s ) 

£ is a K—PathSet 

2.2. Simple Factorization 

Definition 2.5. Consider an edge e = ({v , w}, n) e E of a graph G = (V, E) 
and define the following equivalence relation in V: a ~ b if a = b or {a, b} = 
{v,w}. Consider the suryective canonical function 7r : V — > V/ ~ such that 
7r(a) = [a]^. We define the contraction of an edge e in G as the graph G ■ e 
such that 

G ■ e = {V/ ~, E ■ e) 

where E ■ e = { ({ir(a), 7r(6)}, n) / ({a, 6},n) G £ - {e} } (see Figure |2|. We 
will denote by K e = tt(K) the new distinguished set of nodes contained in 
V(G ■ e) where K is the distinguished set of nodes in V{G). 

Definition 2.6. Consider an edge e = ({v,w},n) G E of the graph G = 
(V, E). We define the deletion of the edge e of G as the graph (see Figure [T]) 

G-e=(V,E-{e}) 

The following is the simple factorization proposed for the first time by 
Moskovitz [Mo]: 

Proposition 2.1. Consider a stochastic graph G and a subset K of its nodes. 
For any edge of G we have that 

R K {G) = PiR Kl {G ■ ti) + (l- P i)R K {G - e<) 
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Figure 1: Deletion of the edge e 




Figure 2: Contraction of the edge e 



Definition 2.7. Define the following partial order in the set of states of G: 
£ < J 7 if C The state £ is a i^-minpath if it is minimal in 

the set of .fT-PathSets. 

The if-PathSets have the following coherence property: If £ < J 7 and £ 
is a .fT-PathSet, then J 7 is a if-PathSet. For every f^-PathSet J 7 there is a K- 
minpath £ such that £ < J 7 . We conclude that the set of i^-PathSets equals 
the set of states J 7 greater or equal to some i^-minpath. This motivates the 
following definition [Co]: 



Definition 2.8. An edge e of a graph G will be called i^-irrelevant if £ (e) 
for each i^-minpath £. 







Proposition 2.2. If e is a K-irrelevant edge of a stochastic graph G, then 

R K (G) = R K (G - e) 

2.3. K- Reliability Polynomial 

Consider a graph G and the stochastic graph G p whose underlying graph 
is G and its Bernoulli variables are identical and independents with parameter 
p. By proposition |2.1| and the finiteness of the graph, we have that the K- 



Reliability is a polynomial respect to the parameter p: 
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Definition 2.9. Rk{p) = Rk{G p ) is the es K-Reliability polynomial of G. 

See that 

n 

r k {p) = Y,c«p\i-pT- 1 

i=0 

where Cf is the number of K-pathsets with exactly i operative edges and n is 
the number of edges in G. In particular, the calculation of the K-Reliability 
polynomial is equivalent to the calculation of the numbers . Each K- 
minpath is a Steiner tree of G which covers K (if K is the set of nodes of G, 
then a i^-minpath is a spanning tree). The calculation of the K-Reliability 
polynomial solves the following iVP-complete problem: Given a number b 
and a subset of nodes K, Is there a Steiner tree which covers K with edge 
number less than or equal to bl To answer this question we just have to see 
weather or not m < b where Cf is the first non zero coefficient of the above 
expression. In fact, we have the following [Co]: 

Proposition 2.3. The Calculation of the K-Reliability polynomial is a NP- 
hard problem. 

3. Combinatorics of the Problem 

Hypothesis 1: In the whole paper, G±, G2 and G are graphs with 
distinguished subset of nodes K\, K 2 and K respectively such that 
G = G x U G 2 , K = K x U K 2 and 

{k u k 2 , • • • k n } = K x n K 2 = Gi n G 2 

Hypothesis 2: We assume in the whole paper that for each node v G K 
there exists a path in G that joins v with some vertex ki G K\ D K 2 . 

These Hypothesis are illustrated in figure [3j In view of the first hypoth- 
esis, it is reasonable to assume the second one, otherwise Rk{G) = and 
there would be no necessity for any calculation. For notational convenience, 
the K subscript in R K (G) will be omitted in the rest of the paper. 

Lemma 3.1. G is K-connected if and only if G is {ki, k 2 , . . . k n } -connected. 

Proof: The direct is trivial. Conversely, take a pair o nodes a and b in 
K. There are paths P a y P& connecting a and b with fc, y kj respectively. 
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Figure 3: Hypothesis on the Graph 



Because G is {k\, k 2 , . . . fc n }-connected, there is a path P connecting fc, with 
kj. The concatenation of the paths P a , P and joins a with b. □ 
The previous lemma motivates the following definition. 

Definition 3.1. Consider the equivalence relation: ki ~* kj if there is a path 
in G\ joining ki with kj, I = 1,2. The connectivity state of Gi is the partition 
of {kx, k 2 , ■ ■ ■ k n } given by 

Denote by Con the set of partitions of {ki, k 2 , . . . k n }. Figure [i] shows 
some useful notational and diagrammatical ways to represent a connectivity 
state. 

Definition 3.2. For each connectivity state A denote by Gf the graph 
resulting of the identification of the nodes in {ki, k 2 , ■ ■ ■ k n } of Gi by the 
state A; i.e. given the graph Gi = (V, E) define Gf = (V/ E A ) where 

E A = { ({7r(o),7r(6)},n) / ({a,b},n) G E } 

and ~^ is the equivalence relation in V generated by A with the canonical 
suryection 

7T ; V -> V/ ~f 

such that 7r(a) = [a]^. The distinguished set of nodes is K A = tt{K{). 

Definition 3.3. The set of connectivity states of G is Con x Con and (Ci, C 2 ) 
is the connectivity state of G where C\ is the connectivity state of Gi, I = 1,2. 

The rest of the section is devoted to the translation of the probabilistic 
problem given in the introduction into a purely combinatorial one. 
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{{1,3},{2,4» = 




Figure 4: Connectivity State 

Definition 3.4. Consider a stochastic graph G = (V,E, (p e ) eg £) and define 

Pi{C) = P(C = Ci) 
P G (A,B) = P((A,B) = {C X ,C 2 )) 

where A, B and C are connectivity states and I = 1, 2. 

The connectivity states can be seen as a partition of the sample space of 
states of Gi where I = 1, 2, so 

i = E p ^ 

CeCon 

Because the random Bernoulli variables of the edges of G are independent, 
connectivity states of G\ are independent of those of Gi so Con x Con is a 
partition of the sample space of states of G and 

P G (A,B) = P 1 (A)P 2 (B) 

Definition 3.5. We say that a connectivity state (A, B) of G is connected 
if 

{k u k 2 ,... k n }/ ~= {{&i, k 2 ,... k n }} 

where ~ is the following equivalence relation in {k\, k 2 , . . . k n }: k t ~ kj if 
ki ~^ kj or ki ~ B kj. 

Figures [5] and [6] show examples of connectivity states of G. This way we 
have the following formula for the i^-reliability of G: 

R(G)= £ E PM)P2{B) 

(A,B) connected (-^iS) connected 
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Figure 6: Non connected connectivity states of G 

also given in [Ro] (under a different notational scheme) where an algorithm 
for the reliability exact calculation is developed based on it. The above 
formula is not a factorization theorem. This formula suggests the following 
algebraic construction: Consider the Q- vector space V generated by Con 
and the linear functional Pi : V — > 1R which is the linear extension of the 
probability Pi (we are abusing of the notation and denoting the functional 
by the same name). This way we get the functional 

Pi (8) P 2 : V <g> V ->• R 

and the following expression for the reliability: 

R(G) = P 1 ®P 2 l £ A ® B\ 

\ {A,B) connected I 

Definition 3.6. Denote by R^(Gi) the i^-reliability of Gf 
Lemma 3.2. R A (G t ) = £ B / (AB) connected Pi(B) 

Proof: Observe that Gf is /^-connected if and only if Q verifies that 
(A,Ci) is connected so 

R A (Gi) = P(Ci e{B / (A,B) connected}) 



8 





+ c + 



Figure 7 



and because the set of connectivity states is a partition of the sample space 
of the states of Gi, we get the result. □ 
We can rewrite the above formula with the functional Pf. 

Ra(Gi) = Pi[ Yl B 

. B I (A,B) connected 

Then we have finally translated the original probabilistic problem in the 
following combinatorial one: 

Problem (Combinatorial Version): Express the vector R 6 V <g> V 
given by 

R= Y A®B 

(A,B) connected 

in terms of the vectors R^ E V given by (see Figure^ 

ra= b 

B I (A,B) connected 

The relation between the combinatorial and the probabilistic problem is 
given by the functional P\®P2 '■ V®V — > R; i.e. Acting on the combinatorial 
expression by Pi <8> Pi gives the probabilistic one. 



4. The Factorization Theorem 

Considering an ordering in the base Con we can write 

m 

R = aijAi <g) Aj 

where m is the cardinal of Con and A = (a^) is the matrix given by = 1 
if (Ai,Aj) is connected and a^- = if it is not. The matrix A will be called 
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the connectivity matrix. Consider the linear operator T : V — > V such that 
T(A) = R A . 

Lemma 4.1. The connectivity matrix is symmetric and is the associated 
matrix of the operator T relative to the ordered base Con. 

Proof: By definition, A es symmetric. Relative to the base Con we have 

[T(Aj)]con = [RAj]con = {o>lj, a 2j, ■ ■ ■ d m j) 

SO 

PI Con = ( a ij) 

□ 

The next proposition will be proved in the next section (Proposition 5.7). 
Proposition 4.2. The operator T has determinant 

det{T) = ± Yl {m A -l)\ 

AeCon 

where m A is the number of classes in the connectivity state A. In particular, 
T is an automorphism ofV; i.e. T G Autq(V). 

Observation 4.1. The operator %T : Z(Con) — > Z(Con) such that T(A) = R_a 
whose relation with the operator T is 

T = Q ® zT 

is not an automorphism of the Z- vector space generated by Con if n > 2. In 
fact, its cokernel is (Proposition 5.10) 

Z(Co? 

Con) A£Con 

where m A is the number of classes in the connectivity state A. 

Observation 4.2. If we write the above proposition in terms of the connec- 
tivity matrix, 

det(A) = ± Y[ (m A -l)\ 

AeCon 

we see that the left hand side of the equality is in some sense topological (it 
is related to connectedness) while the right hand side is combinatorial. It is 
a beautiful result. 
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Finally we have the factorization theorem which solves the posed combi- 
natorial problem. This is the main theorem of the paper. 

Theorem 4.3. Let (6^) = A^ 1 where A is the connectivity matrix. Then 

rn 

and the above expression doesn't depend on the order of the base Con. 



Proof: By lemma 4.1 we have 

so 

m 

R = ^2 a ijbkibhjRA k ® RA h 
h,i,j,k=l 

Because (ay) and (fry) are inverses 

m 
i=l 

and we get 

m m 

R = o~kjbhjRA k ® = bhkRA k ® -SU,, 

h,j,k=l h,k=l 

Finally, because ^4 is symmetric, its inverse is too so 

m 

We have to show that the expression above is independent of the order in 
the base Con. It is enough to make an intrinsic formulation not depending 
on any order of the base: 



R = ^2 A®R A 

A&Con 

= J2 T-\T{A))®R A 

A&Con 

= J2 t ~\ r a)®Ra 



A&Con 
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Figure 8 

The above expression is the one we were looking for. Fixing an order in 
the base, the last expression reduces to the original one relative to this order. 
□ 

As it was mentioned in the previous section, we solved the posed proba- 
bilistic problem just applying the functional P\ <g> P2 on the last expression: 

Corollary 4.4. Under the hypothesis and notation of the previous theorem, 



The case n = 1 is clear and reproduces the well known factorization 
respect to an articulation point. Let's see how the theorem works in a cutset 
of size two; i.e n = 2. Ordering the base Con by 

Con = {12,^2^} 
we get the connectivity matrix 



.4 



and its inverse 



.4 



-1 



1 

1 1 

-1 1 
1 



so 



R = R 



12 



Rfa — R12 ® R12 



R12 + R12 

Figure [8] shows this factorization. See the appendix for another proof of this 
result. 

Let's see the case n — 3. Ordering the base Con by: 



Con = {123,1 23, 13 2, 12 3,123} 
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we get the connectivity matrix 

/ 1 \ 

111 

10 11 

110 1 

\ 1 1 1 1 1 / 

/ 1 -1 -1 -1 2 \ 
-1-11 10 
-11-110 
-11 1-10 
\ 2 / 

The following expression for the n = 3 factorization is illustrated in Figure 

ED 

# = #123 ® #123 + #123 ® #123 + ■ ■ ■ 
• • • + 2^123 ® #123 - #123 ® #123 ~~ #132 ® #132 ~ #123 ® #123) 

+ 2(^123®#132+#123®#123+#132®#123 _ #123®#123 _ #123®#132 _ #123®#123) 
+ 2(^132®#123+#123®#123+#123®#132 _ #123®#123-#132®#123-#123®#123) 

5. The Connectivity Matrix 

5.1. Connectivity Matrix Determinant 

In this section we will identify a connectivity state with its symmetry: 
i.e. there is a one to one correspondence 

S A ->A 

between the connectivity states and the subgroups generated by transposi- 
tions of the permutation group S n such that 

A = {h, k 2 ,... k n }/S A 

This means that the connectivity state A is the set of orbits of {ki, k 2 , . . . k n } 
under the action of the subgroup generated by transpositions S A . 



A 



and its inverse 

2 
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Figure 9 



Observation 5.1. The subgroups must be generated by transpositions to get 
a maximal subgroup under the above condition and to get the inverse corre- 
spondence 

A^S A 

We could think of the connectivity states as representations of the subgroups 
generated by transpositions of the permutation group. 

This identification shows that the connectivity states have a natural com- 
mutative monoid structure with unit 

e = {id}^{{h},{k 2 },...{k n }} 

under the product 

A ■ B = (A U B) 

Because the product of subgroups generated by transpositions y generated by 
transpositions, this product is well defined. Observe that under the identifi- 
cation, Con has also a natural partial order given by inclusion. From now on 
we will denote by the same name Con the set of partitions of {ki, k 2 , . . . k n } 
and the set of subgroups generated by transpositions of the permutation 
group S n . 
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This identification provides a useful and elegant way to characterize the 
connectivity state of G: (A, B) is connected if and only if 

A ■ B = S n 

Definition 5.1. The permutation group S n acts by conjugation in the con- 
nectivity states: 

a (A) = crAa' 1 

where A G Con. 

In particular we have conjugation classes in Con; i.e. the orbits in Con 
under the action of S n : A ~ B if there is a permutation a such that A = 
a(B) = aBa- 1 . 

Lemma 5.1. Consider a permutation a in S n . Then a : Con — > Con is 
unital monoid morphism; i.e. o~({id}) = {id} and 

a(A ■ B) = a{A) ■ a{B) 

Proof: It is clear that a({id}) = {id}. Let's show that a is a monoid 
morphism. See that a(A ■ B) contains a (A) and cr(B) so, by definition, 

a(A ■ B) D a (A) ■ a(B) 

See also that, a^ l {a{A) • cr(B)) contains A y a B so, by definition, 

A-Bca- 1 (a(A)-a(B)) 

Acting by a we get 

a (A ■ B) C a (A) ■ cr{B) 

□ 

This way, the Q-vector space V generated by the commutative monoid 
Con is an associative and commutative Q- algebra with unit. The unital 
monoid morphism a extends linearly to a unital algebra morphism where a 
is a permutation. 

Consider the linear operator n : V — > V such that ir(S n ) = S n and for 
every connectivity state A distinct from S n , 

Tr(A)= JJ (A-A-B) 
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Lemma 5.2. Consider a connectivity state A. The vector tt(A) verifies the 
following properties: 

1. B ■ ir(A) = Vi3 / B it A 

2. C ■ 7r{A) = tt(A) VC / C<A 

3. n(a(A)) = v{ir{A)) Va G S n 

In particular, A ■ 7r(A) = n(A) and B ■ n(A) = for every connectivity 
state B distinct and conjugated to A. 

Proof: 

1. The algebra is commutative and Con is a base of idempotents; i.e. 
A 2 = A for every connectivity state A. 

2. It is enough to see that C • A = A if C < A. 
3. 

a (it (A)) = o-i Yl {A-A-B) 

= II HA) - a(A) ■ a(B)) 

B I BftA 

n (*c*) - ^ ■ ^)) 

<t(B) I trQS)jta(A) 



where we used in the las identity that a is invertible. 



□ 



Definition 5.2. For each connectivity state A we define its connectivity 
number as the coefficient of S n in the expression of 7t(A); i.e. 

tt(A) = . . . + a A S n 



See that, by the third item of lema 5.2, the connectivity number is in- 
variant under conjugation: 

c*A = oc a (A) 

for every permutation a. The object of the following lemmas is to calculate 
the connectivity numbers and show that they are non zero. 
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Definition 5.3. Denote by T n the graph with n nodes and one edge joining 
every pair of nodes. Denote by the resulting graph from the identification 
of the nodes {1, 2, ... n} in T n by the classes of A. 

Let G be a stochastic graph and consider its reliability polynomial R(G). 
Denote by mgr(R(G)) the term of R(G) whose degree equals the edge number 
of G. See that if G has an irrelevant edge, then mgr(R(G)) = and in case 
mgr(R(G)) is non zero, then this term equals the highest degree term of the 
polynomial. 

Lemma 5.3. Consider a stochastic graph G with k edges between a pair of 
distinct nodes i and j of G. Consider the resulting graph G by deleting k — 1 
edges between the nodes i and j ofG. Then, 

mgr(R(G)) = (-p) fc_1 mgr 

Proof: The result is clear for k = 1. Suppose there are k > 1 edges 
between the nodes % and j of G and that the result holds for an amount less 
than or equal to k — 1 of them. Consider an edge a between the nodes i and 
j. A simple factorization on the edge a gives 

R(G) = p R(G ■ a) + (1 - p)R(G - a) 

where G ■ a is the resulting graph by the contraction of a and G — a is the 
resulting graph by deleting a. See that the edge number of G ■ a and G — a 
is the edge number of G minus one and because k > 1, G ■ a has irrelevant 
edges. This way, 

mgr(R(G)) = (— p) mgr(R(G — a)) 
By the inductive hypothesis, we get the result. □ 
Lemma 5.4. The reliability polynomial ofT n is 

R{V n ) = ±(n - 1)! p 9 + . . . 

where the highest degree of the expression g equals the edge number ofY n : 
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Figure 10: Relation between the graphs T 



Proof: In this proof we will make an abuse of notation identifying the 
reliability polynomial with its graph. We claim that 



mgr(T n+1 ) = (-l) n+ n mgr(T n ) p n 

Because T 2 = p and Ti — 1 we have the result for the n - 
the claim is true for every natural less than or equal to n. 



1 case. Suppose 



The Figure 10 shows the relation between the distinct graphs T. By a 
simple factorization on the edge joining the nodes n and n + 1 of the graph 
r, t+ i and the above lemma we have that 



n+1 



pi-pT- 1 T n + ... + (l-p) H, 



n+1 



where the dots denote terms whose degree is less than the edge number of 
T n+ i and the graph H n+ \ results from deleting the edge joining the nodes n 
and n + 1 of the graph r n+ i, see Figure 11 
By the inductive hypothesis, 



n-l 



mgr(T n ) = {-l) n {n - 1) mgr^T^) p 
and the fact that the relation between the graphs T n and r n _! is the same as 



the one between H n+ i and T n (see Figures 10 and 11), we have the following 
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Figure 11: The graph H n+ \ 

relation 

mgr{H n+1 ) = (-1)> - 1) mgr(T n ) p n ~ l 
Then, we have that 

rw = Pi-vT- 1 r„+. . .+(i- P )(-i)>-i) r riP n - 1 +. . . = (-i) n+1 n r nP n +. . . 

where the dots denote terms whose degree is less than the edge number of 
Tn+i- We conclude that 

mgr(T n+1 ) = (-l) n+1 n mgr(T n ) p n 

which proves the claim. This recursive relation shows that mgr{T n ) is non 
zero so it equals the highest degree term of Y n : 

T n = (_i)«+(«-i)+-2(„ - i)t p ("-i)+-i + . . . 

and this concludes the lemma. □ 

Lemma 5.5. Consider a connectivity state A with m classes a%, a,2, ■ ■ ■ a m . 
Then 

R(V£) = ±(m - 1)! p 9 + . . . 

where the highest degree g of the expression equals the non irrelevant edge 
number of T^: 

g = ^2 aia i 

Proof: We will make the same abuse we did in the proof before identifying 
the reliability polynomial with its graph. The graph has m nodes (these 
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are the m classes of A), y ^ J irrelevant edges in each node i respectively 

and diOj edges joining the nodes i and j. 

Consider the graph f ^ resulting from deleting all the irrelevant edges of 
the graph T^. This way and have the same reliability polynomial. 



By the lemma 5.3 we have the following relation between the graphs and 



mgr(T^) = (-p)^^- 1 ) mgr(T m ) 

= (-p)E^-(oi^-i) ( ± ( m _ p 



m 
2 



and this concludes the proof. □ 
Lemma 5.6. Consider a connectivity state A with m classes. Then, 

cla = ±(m — 1)! 

Proof: We claim that 

n(A)=A-l n ( e "^» 

\t transp. / (r)jiA 

In effect, consider a connectivity state B such that B A. There is a 
transposition r in £> not belonging to A. Because A C A ■ B and (r) C A - B 
we have that ^4 • (t) C ^4. • B so 

.A • B ■ (A - A • <r» = A • B - A ■ B = 

This implies that 

AB-( n (A-A-(r))\=0 

\t transp. / (r)-ftA I 

so we have the following expression for the vector ir(A): 

n(A)= J] (A-A-(r})=A-l J] (e-(r» 

t transp. / (t)^.A \t transp. / {r)ftA. 
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which proves the claim. In particular, the last expression implies that 

ola = Co - Ci + C 2 - C 3 + . . . 

where Cj is the number of subsets F with cardinal % of the set of transpositions 
{r transp. / (r) <jt A] such that (F) ■ A = S n . 

Identifying the transposition (ij) with the edge joining the nodes i and j 
of r n , it is clear that Ci is the pathsets number (operational states number) 
of with just i operational edges. This way we have that 

#(0 = c (i- P y + c 1 p(i-py- 1 + c 2 p 2 (i- P y- 2 + ... 

= (-iy(C -C 1 + C 2 -C 3 + ...)p 9 
= {-l) 9 a A p 9 + ... 

where (7 is the non irrelevant edge number of and the dots denote terms 
whose degree is less than g. By the above lemma this implies 

a A = ±(m - 1)! 

□ 

Consider the connectivity state conjugation classes {Oi,0 2 , . . .Ok}. The 
number of classes a connectivity state has is invariant under conjugation so 
we can define mo, as the number of classes in some connectivity state in 0,. 

Identifying each level of the Hasse diagram (inclusion diagram) of the 
subgroups generated by transpositions of S n and because conjugated states 
necessary belong to the same level of this diagram, we can order the connec- 
tivity states in Con in the following way: We order some conjugation class Oi 
of the first level, then we order some other conjugation class Oj of the same 
level and we continue the process until we have order all the conjugation 
classes of the first level. After that, we order the second level in the same 
way as we did in the first and so on until we have order all the connectivity 
states. The given order in Con verifies 

Ai< Aj i < j 

We will call a coherent order an order in Con verifying the above property. 

Proposition 5.7. The determinant of the operator T such that T(A) = R A 

is 

k 

de*(T) = ±J](m 0i -l)! i0i 

i=l 
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Proof: Choose a coherent order in the base Con and consider the con- 
nectivity matrix A relative to this order. Consider the matrix B associated 
to the operator tc relative to the chosen coherent order in Con. Because of 
the identity proved in the last lemma 

n(A)=A-\ J] (e-(r» 

t transp. / (r)ytA- 

we have that B is an inferior triangular matrix with ones in its diagonal, 



B 



( 1 ... \ 
* 1 ... 



V 



In particular, 



det(B) = 1 

We can think about the expression of the vector tt(A) in terms of the base 



Con as elementary row operations on the matrix A so, by lemma 5.2 
have the following identity: 



we 



B l A 



* ae>2 ho-i 



\ 



o \ 



• OL 0k I i0k J 



Because de^B 1 ) = det(B) = 1 we conclude that 



det(A) = Y\ a Oi 



to, 



i=i 



and by the last lemma we finish the proof. □ 
As an example, let's see the n = 4 case. The connectivity states conju- 
gation classes are: 



22 



o 5 



o 2 



o 3 



o 4 




then we have that md equals 4, 3, 2, 2, 1 respectively. By the above proposi- 
tion the determinant of the connectivity matrix is: 



5.2. Connectivity Matrix Inverse 

We can give another monoid structure on Con, defining the product as 
the intersection of the subgroups generated by transpositions where S n is 
now the unit. In the same way as before, each permutation in S n is a unital 
monoid morphism and its extension to a the Q-algebra V (commutative 
and associative with unit, generated by the monoid Con) is a unital algebra 
morphism. 

Consider the linear operator £ : V — > V such that £( {id} ) = {id} and 
for every connectivity state A distinct from {id} we have 



Lemma 5.8. Consider a connectivity state A. The vector ^(A) verifies the 
following properties: 



det(A) = ±(4-1)! 



(3-1)! 6 (2- 1)! 3 (2- 1)! 4 



(1-1)! 



±384 



C I C<A 



1. Bn£{A) = WB I B < A 

2. Cn £(A) = £(A) VC / C>A 

3. £(a(A)) = a(£(A)) W e S n 



Proof: The same proof as in lemma 5.2 



□ 
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Theorem 5.9. Consider a coherent order in Con and the matrices B and 
D associated to the operators tc and £ in the base Con. Consider also the 
matrix 



( cr \ i i0l o . . . o \ 



c 



a ol hoi 



Then, 







A' 1 = B C D 







Proof: Following the proof in 5.7, by lemma 5.2 ij component of the 



matrix C B A is one if Aj < A%\ i.e. Aj r\A% = A% and zero otherwise. 

We can see the expression of the vector £(*4) in terms of the base Con 
as elementary row operations on the matrix C B l A and observing that this 



expression only has elements C such that C < A, by lemma 5.8 and the fact 
that B A and C < A implies B <jt C, we have that 

I = D* C B l A 

Because A is symmetric, the transpose of the last expression shows that A 
has also a right inverse so A^ 1 = B C D. □ 
As an example consider the n = 3 case. A coherent order on the base 
Con is 

Con = {123, 1^,^2,^3,^23'} 
Let's Calculate the matrices B, C y D of the previous theorem: 



tt( 123 ) 
?r( ^3 ) 
?r( ^T2 ) 
tt( ll£T ) 
tt( ^2? ) 



123 - 12 3 
'1^3 - ^23 
^13^2 - ^23 
1^23^ - / 12? 



13 2 



1 23 



2. 123 



24 



and the matrix associated to n on the base Con is 



B = 



( 1 \ 
-1 1 

-1 1 

-10 1 

\ 2 -1-1-11/ 



See that the entries of the last row are the respective connectivity numbers 
so 



C 



( \ 



\ 



-1 



\ 



1 / 



a 123 ) 

£(^3) 
£('13*2) 
C(l"23") 
£( ^23) 



123 

'13*2 
l^T 

m - 



123 
123 
123 
"12*3 



13 2 - 1 23 + 2. 123 



and the matrix associated to £ on the base Con is 



D = 



(l -1 
1 



V 



-1-1 2 \ 

0-1 

1 -1 
1 -1 

1 / 



By the previous theorem, the connectivity matrix inverse is 



/ 1 



A- x = 



1 





1 





\ 



1 / 



(\ 



\ 



1 / 



/ 1 



-1 -1 



1 

1 



2 \ 

-1 
-1 
-1 

1 / 
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The permutation group S3 has the peculiarity that its minimal subgroups 
coincide with its maximal ones and this is the reason why the matrix D 
equals the transpose of B in this particular case. This phenomenon occurs 
only in this case. The inverse connectivity matrix calculation follows exactly 
the same steps as in the case before and it is realized through the algebra 
generated by the connectivity states monoid. 

5.3. Invariant Factors 



By the theorem 5J3 we can see the connectivity numbers unless sign as 
the invariant factors of the abelian torsion group: 

Z(Con) 
Z(R A I A G Con) 

In effect, we have the following theorem: 

Theorem 5.10. 

k 

9 Z K-1)! 



Z(Con) 



Z(R A I A G Con) 



Proof: Following the proof of theorem 5J3 because B and D are trian- 
gular matrices with integer entries and ones in its diagonal, B and D are 
invertible in M m (Z); i.e. 

B,D G GL m (Z) 
In particular £ is an automorphism of Z(Con): 

£ G Aut(Z(Con)) 

This way we have that 

/ a 0l I i0l ... \ 
a 02 I$o 2 ■ ■ ■ 



DA 



B- 1 



\ ... a 0k ho k ) 

Consider the subgroup iV of Z(Con) generated by the vectors a A A such 
that A is a connectivity state. Because of lemma |4.1 , the fact that B is 
invertible in M m (Z) and the identity 

m m m m 

£,{Ra 3 ) = ^aij^i-Ai) = ^2 d ki(iijA k = 22 a k^kiGijA k = 22e ij (a i Ai) 



i,k=l 



i,k=l 



i=l 
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where (e^) = B^ 1 and (dy) = D, we conclude that A" is isomorphic by £ to 
the subgroup of Z(Con) generated by the vectors R A . We have the following 
commutative diagram with exact rows: 







Z(R A I A e Con) Z(Con) 

€. - 5. - 
AT > Z(CW) 



Z(Con) 



5(R.A / .Ae Con) 
I 
I 

Y 



i 



=1 ^(">t), 



-D! 







where lemma |5.6| was used. By the five lemma (well known result in Homo- 
logical Algebra, see [We]), we have the result. 
As an example, in the n = 3 case we have 



□ 



Z(Con) 



while in the n 



Z(R A / Ae Con) 
4 case we have 
Z(Con) 



Zo 



Z(R A I A e Con) 



Z fi ©Z° 



Z, 



In contrast to the n = 3, 4 cases, the invariant factors in the n > 4 cases are 
not just the prime decomposition of the connectivity matrix determinant. 
For example, consider the n = 5 case: 






#0 


m 


[12345] 


1 


5 


fl2^345] 


10 


4 


^T^T5] 


15 


3 




10 


3 


^234 5] 


5 


2 




10 


2 


0M5\ 


1 


1 



and we have that 
Z(Con) 



Z(R A / Ae Con) 



Z, 



7 io 



z? 



z, 



1% 



Z: 



35 
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Appendix A. Another Proof of the n=2 case 

In this section, connected means edge-connected. 

Lemma Appendix A.l. Consider a graph G and subgraphs G\ and G2 
with distinguished sets of nodes K\ y K 2 respectively such that G = G\ U G2 
and {a, b} = K\ fl K 2 = Gi fl G 2 (see the next Figure). Define K = KiU K 2 . 
If G is K -connected, then G\ is K\-connected and G 2 is Reconnected, where 
Gi and Ki is the graph and set resulting from the identification of the nodes 
a and b in the graph and set Gi and Ki respectively. 




Proof: Consider k in K\. There is a path in G connecting k with a. 
Then, there is a subpath in G\ connecting k with a or b; i.e, there is a path 
in G\ connecting k with a = b. This is true for every k in K\ so G\ is 
.fTi-connected. In the same way as before, G 2 is ^-connected. □ 

Lemma Appendix A. 2. Under the hypothesis of the previous lemma, G 
es K-connected if and only if one of the following items holds: 

i. G\ is K\-connected and G 2 is K 2 -connected. 

ii. G 2 is K 2 -connected and G\ is K\-connected. 

Proof: By the previous lemma G\ is i^i-connected and G 2 is ^-connected. 
Suppose that G 2 is not ^-connected; i.e., there are nodes u, v in K 2 such 
that there is no path in G 2 connecting them. 

We claim that Gi is Ki-connected. Suppose that it is not; i.e., there are 
nodes x, y in K\ such that there is no path in G\ connecting them. However, 
G is K-connected through paths 71 and 72 in G connecting x with u and y 
with v respectively. Because 71 and 72 contain the nodes a or b (but not the 
same node), we conclude that a or b are not connected in G. In effect, if they 
were connected by a path in G, then it would exist a subpath connecting 
them in G\ or G 2 so x and y or u and v would be connected. This is absurd 
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because G is if-connected and the nodes a and b are in K. We conclude that 
G\ is A^-connected. 



Conversely, suppose without loss of generality that G\ is fTi-connected and 
G 2 is A 2 -connected. We claim that G is A-connected. In effect, consider ki 
and k 2 in K. If k\ and k 2 are in Ki, then there is a path in G\ connecting 
k\ with &2 and because G\ G G this path is also in G. If k\ is in Ai and 
k 2 is in A2, then there is a path in G 2 connecting k 2 with a or b. Suppose 
without loss of generality that k 2 is connected with a. There is a path in G± 
connecting a with k\. The concatenation of these paths connects k\ with k 2 
in G. The argument for the case k\ in K 2 and &2 in Ki is the same. Finally, 
if fci and k 2 are in K 2 , then there are paths in G 2 connecting k\ and k 2 with 
a or b. If these paths connect k\ and &2 with the same point, then k\ and k 2 
are connected in G2 C G, otherwise there is a path in G\ connecting a with 
b and the concatenation of these three paths connects k\ with k 2 in G. □ 
As a corollary we have the factorization formula in the n = 2 case (see 
figure [8]) . 



Theorem Appendix A. 3. Under the hypothesis of lemma Appendix A.l 
we have that 

Rk(G) = Rk^G^R^G^ + R^G^Rk^G^ — R Kl (G 1 )R K2 (G 2 ) 



Proof: By the previous lemma |Appendix A.2[ the set of i^-PathSets is 



the union of the sets of states C\ y C 2 which verify the items i and ii of the 
lemma. The intersection of these sets is the set of states verifying that G\ is 
Ki-connected and G 2 is ^-connected. From the identity 

P(A UB) = P{A) + P(B) - P(A n B) 

and the independence follows the result. □ 
As a corollary, taking a = b in the previous formula we get the well known 
factorization respect to an articulation node. 
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Appendix B. The n=4 case 

For the case n = 4, ordering the base Con by: 



Con = {1234, 12 34 ,13 24 , 23 14 , 1 23 4, 13 24, 12 34, 14 23 , 
. . .^^^4^,^2^^3^,1^3?,2^3?,31^?,1^3 S 4,T23^} 

we get the connectivity matrix 













































1 \ 























1 


1 











1 


1 


1 























1 





1 





1 





1 


1 


























1 


1 


1 








1 


1 


























1 


1 





1 


1 





1 























1 





1 


1 





1 





1 























1 


1 





1 


1 








1 





1 


1 








1 


1 





1 


1 


1 


1 


1 


1 


1 





1 





1 


1 





1 


1 





1 


1 


1 


1 


1 


1 








1 


1 


1 


1 





1 


1 





1 


1 


1 


1 


1 











1 





1 


1 


1 


1 


1 





1 


1 


1 


1 








1 





1 





1 


1 


1 


1 


1 





1 


1 


1 





1 








1 


1 





1 


1 


1 


1 


1 





1 


1 





1 


1 


1 











1 


1 


1 


1 


1 


1 





1 


1 
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1 


1 


1 
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1 


1 


1 


1 


1 


1 


1 


1 


1 / 



and its inverse 
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/ 


_1 

X 


i 

X 


i 

X 


i 

X 


i 

X 


1 

X 


1 


_i 

X 


_1 

X 


_1 

X 


_2 

Zi 


_2 

Zi 


_2 

Zj 


_2 

Zi 


6 \ 

U i 




1 

X 


2 

Zi 


_i 

X 


_i 


_L 


_L 




1 

X 


1 

X 


_2 

Zi 


_1 

X 


_1 

X 


2 

Zi 


2 

Zj 


n 




1 

X 


_1 

X 


9 

Zi 


_1 

_L 


_1 

X 


_1 


X 


1 

X 


_2 

Zi 


1 

X 


X 


9 

Zj 


_1 

X 


2 

Zi 


n 




1 

X 


_1 

X 


_1 

X 


2 

Zi 


_1 

X 


_1 

± 


X 


_2 


1 

X 


1 

X 


2 

Zj 


_1 

X 


X 


2 

Zj 


n 




1 


-1 


-1 


-1 


2 


-1 


-1 


-2 


1 


1 


-1 


2 


2 


-1 







1 


-1 


-1 


-1 


-1 


2 


-1 


1 


-2 


1 


2 


-1 
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-1 
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-1 


-1 
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-2 
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-1 


-1 


-1 


1 
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1 


-1 


-1 


-1 


-1 





V 


6 












































^ — -^1234 ® "^1234 + -^1234 <S> -R1234 + • • • 
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Appendix C. Factorization of the Random Cluster Model in the 
non Cluster Limit 



In Statistical Mechanics, the random cluster model [Gri] partition func- 
tion of a stochastic edge connected graph G = (V, E, (p e ) eeE ) is given by 

Z & G ) = E II V kH Pe {e \l - Pe) 1 ^ 
uiGCl e£E 

where Q = {0, 1} E is the set of states and k(u) is the number of connected 
components of the state u. The parameter q gives weight based on the num- 
ber of closed sets of connected vertices including isolated ones; i.e. clusters. 
The q = 1 case is the percolation model and the q > 1 case prefers more 
clusters while the q < 1 case prefers fewer clusters. The non cluster limit is 
the q ~ case. It is easy to see that a factorization of G as the one described 
in the paper is realized as the product of partition functions in percolation 
theory: 

Z(1,G) = Z(1,G 1 ) Z(1,G 2 ) 

and this is the only known factorization in the cluster random model. Be- 
cause of the following relation between the random cluster model and the all 
terminal reliability: 

Z(q,G) = qR All (G) + 2 (q) 



we have the following factorization in the non cluster limit (see Theorem 4.3 ): 



dZ(q,G) l A 9%Gf) dZ(q,G^) l 

dq U =° ~ 2-, ^ dq l*=° dq ^ 



where b is the connectivity matrix inverse. 
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